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Abstract 

Background: Based on a literature review, the current study aimed to construct 
mathematical models of lactate production and removal in both muscles and blood 
during steady state and at varying intensities during whole-body exercise. In order to 
experimentally test the models in dynamic situations, a cross-country skier performed 
laboratory tests while treadmill roller skiing, from where work rate, aerobic power 
and blood lactate concentration were measured. A two-compartment simulation 
model for blood lactate production and removal was constructed. 

Results: The simulated and experimental data differed less than 0.5 mmol/L both 
during steady state and varying sub-maximal intensities. However, the simulation 
model for lactate removal after high exercise intensities seems to require further 
examination. 

Conclusions: Overall, the simulation models of lactate production and removal 
provide useful insight into the parameters that affect blood lactate response, and 
specifically how blood lactate concentration during practical training and testing in 
dynamical situations should be interpreted. 

Keywords: Aerobic power, Anaerobic power, Blood lactate, Cross-country skiing, 
Muscle lactate, pH 



Background 

The metabolic power in humans is based on the production and consumption of ade- 
nosine triphosphate (ATP). Despite the approximately 100-fold increase in ATP utiliza- 
tion from rest to maximal-intensity exercise, the energetic demands of the muscles are 
usually satisfied without depleting the intracellular ATP e.g., [1-3]. In this connection, 
three sources for ATP synthesis are available; First, ATP can be produced aerobically 
in the mitochondria by oxidative phosphorylation. Second, ATP can be produced by 
anaerobic synthesis due to glycolysis or glycogenolysis. Finally, ATP can be produced 
by phosphocreatine (PCr) break down to Creatine (Cr) (i.e., ADP + PCr gives ATP + 
Cr in the Creatine Kinease (CK) reaction) e.g., [1-3]. 

The rate of oxygen (0 2 ) consumption can be set to the sum of 1) a constant rate 
(resting 0 2 consumption), 2) a rate due to unloaded body movements and 3) a rate 
proportional to the aerobic energy used to perform work. For moderate constant work 
rates, the aerobic power increases towards a steady state condition. The concept of 
maximal lactate steady state (MLSS), that is the highest intensity where a steady state 
lactate can be obtained, has been regarded as important for endurance performance e. 
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g., [4-6]. For exercise intensities above MLSS, associated with sustained acidosis, a slow 
component delays the attainment of a steady state value and causes 0 2 uptake to 
increase to values greater than those predicted from aerobic steady state demands. For 
exercise intensities exceeding the maximal oxygen uptake, the steady state corresponds 
to the level that would be attained if it was possible to carry out the exercise under 
pure aerobic conditions [7]. Obviously, this virtual steady state is never reached as the 
increase in oxygen uptake ends when the maximal oxygen uptake is achieved. 

When the rate of ATP production by oxidative sources becomes insufficient, high 
rates of glycolytic or glycogenolytic ATP production are required. The endpoint of gly- 
colysis is pyruvate, which represents a metabolite that can be reduced to lactate or oxi- 
dized to C0 2 or H 2 0. Thus, by increasing exercise intensities, the working muscles 
and various tissues produce more lactate and release it into the plasma. At the same 
time, the skeletal muscles, the heart, the liver and the kidney cortex remove lactate 
from the circulation, and lactate is suggested to act as an intermediate for the shuttling 
of carbohydrate from cells and tissue with relatively low oxidative capacity to cells and 
tissues with high oxidative capacity [8-11]. Thus, it is well established that the blood 
lactate concentration is the result of the production and the removal of lactate in the 
blood. During steady state sub-maximal exercise, when lactate production (influx) 
equals lactate removal (outflux), the lactate concentration in the lactate pool stays con- 
stant and the rate of oxygen consumption is the measure of the whole body energy 
expenditure regardless of the magnitude of lactate production and removal or the 
absolute blood lactate concentration. At exercise intensities above steady state, a rise in 
the concentration could be attributed to an increase in the rate of lactate production 
or result from a decrease in the rate of lactate removal. 

Lactate itself does not lead to muscle fatigue at high exercise intensities, and fatigue 
is most likely a result of decreased muscle pH and the associated reactions e.g., 
[12,13]. For the skeletal muscles, the pH is normally around 7.1 but can fall to 6.4 dur- 
ing heavy exercise [14], and numerous experiments suggest a negative relationship 
between decreased pH and the muscle contractile function [15-20]. In the capillary 
blood pH is normally (in steady state conditions) around 7.45, but can fall to around 
7.05 during heavy exercise [21]. Previous researchers have found a close relationship 
with a time delay between lactate concentrations in muscle and blood [22]. 

Because muscle groups often work unequally, heterogeneity with regard to lactate 
concentration is likely in muscles during dynamic exercise. A rigorous estimation 
would require application of several compartments with production, removal and 
exchange between compartments. As an approximation, a two-compartment model 
can be applied, where the muscles and other organs that can remove lactate (such as 
the heart, the liver and the kidney cortex) are regarded as one compartment and the 
blood space and other tissues are grouped into a second compartment [23-25]. These 
studies indicate that the flux of lactate into a compartment depends upon the lactate 
gradient and its permeability. Furthermore, blood lactate recovery curves from muscu- 
lar exercise can be described by a bi-exponential time function and a two-compart- 
ment model consisting of the previously worked muscle and the remainder lactate 
space. The time constants of the bi-exponential time function fitted to the arterial 
blood lactate recovery curves reflect the abilities to exchange lactate between the two- 
compartments and the ability to remove lactate from the total lactate space including 
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the working muscle department [26]. The two time constants are found to decrease 
with work rate and duration of the preceding exercise [27,28]. 

The maximum anaerobic energy that can be utilized is proportional to the sum of Cr 
and lactate that can be accumulated in the body. PCr is an energy buffer that supports 
the transient failure of other metabolic pathways to support ATP. The equilibrium 
constant of the CK reaction is around 20 and the slightest drop in ATP allows the 
reaction to proceed to ATP. Thus, the ATP concentration stays nearly constant until 
almost all the PCr is utilized. The PCr levels follow an exponential time course after 
changes in work rate before approaching a steady state condition at moderate exercise 
intensities. In such cases, a strong similarity has been reported for the time constants 
of the 0 2 kinetics and the PCr consumption [29]. However, for exercise intensities 
above lactate threshold, the anaerobic glycolytic energy supply becomes significant and 
the association between PCr and 0 2 rate has not yet been systematically reported. Dur- 
ing recovery, the level of PCr must be recovered, the pH must be re-established and 
ADP removed. While the PCr recovery is mainly due to oxidative ATP synthesis, the 
PCr stores may be rebuilt by anaerobic glycolysis [30-32]. 

The current study aimed to derive mathematical models for production and removal 
of lactate in the blood and muscles during dynamic whole-body exercise. As cross- 
country skiing is a whole body exercise where athletes both train and compete on vary- 
ing terrain and at constantly varying speeds and work rates [33-35] this locomotion 
was used for that purpose. Thus, the mathematical models were compared with experi- 
mental results from the laboratory during treadmill roller skiing. We hypothesized a 
mathematical two-compartment model of lactate production and removal to accurately 
predict blood lactate concentration during steady state and at varying exercise 
intensities. 

Methods 

Overall design 

Initially, the current study derived mathematical simulation models of lactate produc- 
tion and removal in the blood and muscles by utilizing Mathematica 8 (Wolfram 
Research Inc., Champaign, IL, USA). Thereafter, the simulations were compared with 
experimental data from an elite skier performing laboratory tests while treadmill roller 
ski skating both during a steady state and at varying exercise intensities (see details 
below). 

Steady state aerobic power 

aw i 4. aw , i , j .,, . = mvv +iimgCos{a)v + mgSin{a)v u . - 
Work rate W on the treadmill is v J v J , where v is the 

treadmill velocity, (A the coefficient of friction, m the mass of the skier and a ~ Sin(a) 
is the treadmill incline in radians. The power due to the change of kinetic energy mvv 
is zero on the treadmill since the velocity is constant. firngCos(a)v is the power of 
roller friction, and mgSin(a)v is the power of gravity due to the inclination of the 
treadmill. We define Q max to be the maximal aerobic power and Q a to be the steady 

state aerobic power that has been found while treadmill roller skiing. Q a is regarded 
as the virtual steady state aerobic power, with Q a m = d Min (Q a , Q ma x) > where Min is 
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the minimum function. The Min function is employed to ensure that the aerobic 
power does not exceed the maximum aerobic power. Below Q max the virtual steady 
state power Q fl _ Q a is found to be linear with work rate for a given cycle rate and 

incline, and as a hypothesis we apply this linearity also for metabolic powers above 
Q max and find 

ha m = 0^ + Qulif) +C 2 xA(a)x W, Q a m = Min (4, Qmax) 

res t unloaded loaded (1) 

def 

Q b = 80//5, Qulif) ^ UlJ/s,c 2 = 5.8, A(a) = 0.92(1 + 1.19Exp(-71.8a)) 

Qt> is the metabolic power at rest (set to 80 J/s), Q un the metabolic power of 
unloaded movements (zero work rate) which is dependent of the cycle rate. We define 
the cycle rate to be constant in this article. Thus, Q un = 111 J/s here. 

Thereafter, we let Q a be the aerobic power, Q a the virtual aerobic power and 

Q a m = d Min(Q max , Qa{t))' A first order differential equation of aerobic power with the 
virtual aerobic power as input was used to account mathematically for the delay in 
aerobic power with a time lag during steady state work rate. Thus, 

~ A f . mod Q a — Qa(t) f . „ mod . .. , _ Aaaaata 

QflM = — / (a), Qa = Min{Q max ,Qa{t)),{b),Ta= (2) 

305, Q a {t 0 ) = Q b 

The "dot" means time derivative, r is a time parameter quantifying the time before 
the aerobic power reaches steady state during sub-maximal work rates. We uses r = 30 
s according to di Prampero [36]. 

Figure 1 shows the steady state aerobic power as function of the work rate for the 
G3 skiing technique. 

Anaerobic power and lactate concentration in blood and muscle 

Lactate concentration in the lactate pool (C(t)), i.e. the mass of lactate per unit volume 
of this pool, increases only if the rate of lactate appearance (influx) in the lactate pool 
is larger than the rate of lactate disappearance (outflux). The current study uses a 
modified version of Brooks [8] and determines the levels of lactate pool concentration 
by modeling both the influx and outflux streams of lactate and simulate the blood lac- 
tate concentration according to Moxnes and Hausken [37]. Here two pyruvate mole- 
cules are produced for each glucose or glycogen molecule during glycolysis/ 
glycogenolysis. One molecule of pyruvate gives one molecule of lactate. The increase 
in glycogenolysis/glycolysis, due to increased exercise intensity, ends when a maximal 
rate of glycogenolysis/glycolysis is achieved [37]. Therefore, the rate of pyruvate 
appearance (P) has a least upper bound, which is denoted as P max . By neglecting rates 
due to changing pyruvate concentration in the plasma, the rate of lactate appearance 
in mmol/L is R a = P. However, as pyruvate can be oxidized in the mitochondria we set 
the influx of pyruvate into the mitochondria to be a 0 Tanh(P 0) P)/(Po, cc 0 < 1, where a 0 
and p 0 are to parameters fitted to the data and Tanh() is a function that accounts for 
saturation at high lactate concentrations according to Moxnes and Hausken [37]. 
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Figure 1 The aerobic power Q(W, a) as a function of work rate (W) for inclines of a = 0.05 (upper 
line) and a = 0.12 (lower line) during treadmill roller skiing using the skating technique. The curve 
fittings are based on least square fit to the data. Maximal metabolic rate (Q max = 1886 J/s) is 
represented by a straight horizontal line. The lactate threshold (Q LT = 1755 J/s) is represented by the 
straight dotted horizontal line. ■: Experimental values for a = 0.05, A : Experimental values for a = 0.12. 



Thus, the rate of lactate appearance is R a = P - a 0 Tanh(P 0 P)/l3 0 where we expect that 
a 0 is around 1. 

During severe exercise, glycogen re-synthesis by the liver is severely depressed. As a 
hypothesis, we forecast that the rate of lactate disappearance due to both glycogen re- 
synthesis and lactate oxidation is R d = d 0 x (Tanh(x C )lz) x D(Q a ) x (Q max - Q a ), 
where d 0 and % are two parameters that are fitted to the data. D(Q a (t)) is an unknown 
function that is monotonically increasing with the aerobic power. When aerobic power 
equals the maximum aerobic power (i.e. Q a = Q max ) no lactate disappearance takes 
place. Altogether, for a one compartment model of lactate we have 

c(t) = 

Ra ~ Rd m =P - a 0 Tanh{p 0 P/) /#>- 

' * ' (3) 

appemnce v ' 

do x {Tanh{ X C)lx) x D(Q fl ) x (Q max - Q a ) 

' v ' 

disapperance 

In equation (3) we need a model for the rate of pyruvate appearance (P). To our 
knowledge, no such model exists in the literature, and we hypothesize that this func- 
tion is linear with work rate up to P max which is the least upper bound. This simple 
assumption might constitute a potential weakness in our argument. However, if we 
have experimental support this assumption also seems like the most reasonable one. 

Thus, we define a virtual steady state ((P)) analogous to the virtual aerobic steady 
state (Q a ) and hypothesize that 
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p (t) mo d P-m ran = WSiP ™0 Min{Pmaxi p it)) (4) 
Tan 

where z an is the time constant for full activation of glycolysis/glycogenolysis during 
muscle contractions, set to 10 s [38]. P(t) is the steady state rate of pyruvate appear- 
ance. Note that the model in equation (3) applies only for the chosen type of exercise 
and a fixed concentration of glycogen in the body. 

For aerobic power we assume that q a is a linear function of work rate and forecast 

that the steady state rate of pyruvate appearance is approximately proportional to the 
virtual steady state power leading to 

ht) m =P0~Qa (5) 

Due to the rather fast response time in equation (4), an approximation is that 
(t) « P = p 0 Q a (t) • For equation (3) this gives that 

C(t) » po (Q a {t) ~ J^ T «nh (A)poQ«(t))) - d 0 x {Tanh{ X C{t))lx) x D(Q a (t)) x [Q max - Q fl (t)), 
C(t 0 ) = 0.045feg/m 3 = 0.5mmol/L 

Finally, a model for D(Q a ) is needed. As a hypothesis we set that 

D(Q a ) m = oc (q„ - ^-Tanh(p oPo Q a )j (7) 
Here the constant of proportionality can be scaled by d 0 . 

A general solution for constant aerobic power is feasible in (5)-(7) when we let 

Tanh{ X C{i))lx « C(t) (Q fl = Q fl = constant) [37]. 

Altogether, the steady state lactate concentration for constant aerobic power at 
intensities below the MLSS from equation (5) -(7) is 

7 i + x 

/ Qmax) 



C = ArcTanh ( P -^^-) /x = Un + d ^~^ 



Ix 

/Qmax) / (8) 



LimC = Po 



x^o d 0 (Q max -Q a ) 



The maximum aerobic power that can be used for steady state concentration is given 
by 

1 P -^^=0^Q a = Q LT = Q max - P -^^^ = l--^- (9) 

do [Qmax — Qa) ^0 Qmax doQ max 

where Q LT = 0.93Q max = MLSS = , p 0 = 10" 5 'kg/ \m 3 s)l (/ /s), d 0 = 7. 2410" 8 / {J /sf/s, X = 0.95m 3 /kg 
and Q m ^ = 1886 J/s. 

According to equation (8) all steady state lactate concentrations can be achieved for 
steady state aerobic powers below Q LT (MLSS). Furthermore, the steady state lactate 
concentration (C) approaches infinity when the aerobic power approaches Q LT . How- 
ever, it has been discovered that not all steady state lactate levels are tolerated over 
time. This means that blood lactate levels above a certain level of exercise can be 
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terminated before steady state is reached. The anaerobic power due to anaerobic glyco- 
genolysis/glycolysis (Q^ n (0) can be calculated from the increase lactate concentration 
when using the relation from di Prampero and Ferretti [7] as 

QG ( t ) = m x X x C{t), X = 3 x 20 : J/(kgmmol/L) (10) 

Blood lactate concentration usually continues to rise a short period of time after 
exercise. The model in equation (6) does not capture this phenomenon since we only 
considered one compartment; the lactate pool. In the next step we consider the differ- 
ent compartments involved: working muscles, blood and other tissues such as the liver, 
kidney and heart. For exercise powers significantly above resting values, we regard a 
two-compartment model as sufficient: a) the blood compartment and b) muscles and 
other tissues as one compartment (denoted muscles in the rest of current manuscript). 
C b (t) is the lactate concentration in the blood and C m (t) the lactate concentration in 
the muscles. Vb(t) and V m (t) are the volumes of muscles and blood, respectively. The 
total lactate pools volume is V = Vb(t) + V m (t), set to 0.18 L per kg body mass as an 
approximation for the skier modelled here [39]. The muscle mass is set to be 10 kg, 
based on an iDexa scan of the skier, and muscle volume as 10 L. We assume that lac- 
tate moves between the muscle and blood compartment with some time dynamics. 
Thus, we propose the following model for the muscle's (C m ) and blood's (Q>) concen- 
trations of lactate 

CmM = \poD (4(t) - d 0 {Tanh{xC m {t))/x)D{Q a {t)){Q max - Qjl - k la CwW ~ Cb(0 

\V m / L \ / J V m , v 

Cb(t) = kia 77 

I<i a is a parametric function that scales the movement of lactate into and out of the 
blood. The rate of lactate appearance in the blood is set proportional to the difference 
in lactate concentration between blood and muscles. To account for different rates of 
transport into or out of the blood we let K Xa be dependent on C m (t) - C b {t). Estimates 
of these parameters were found by visual curve fitting. Visual curve fitting gives plausi- 
ble values for the parameters based on plotting of experimental data that is compared 
with simulations. The parameters were sought to have biologically trustworthy numeri- 
cal values and a least square fit of the data was performed to produce best fit esti- 
mates. The parameters employed here are 

p 0 = 10- 5 ^/(m 3 5)/(//5), d 0 = 7.2410- 8 /(//s) 2 A, x = 0.95m 3 /kg, a 0 = 0.9 
l/(A>po) = 0.6Q max ,a 0 = 0.9, fen = 0.05/s x L, V b = 4L,V m = 10L 



Experimental tests 

The derived mathematical simulations are compared with experimental data from an 
elite crosscountry skier while roller skiing on a treadmill employing the skating G3 
technique. The mass of the skier was m = 77.5 kg and the friction coefficient on the 
treadmill was ft - 0.024 in all tests. Equipment and procedures were similar to the stu- 
dies by Sandbakk, Holmberg, Leirdal and Ettema [40,41]. All treadmill tests were per- 
formed on a 6 x 3 m motor-driven treadmill (Bonte Technology, Zwolle, The 
Netherlands). Inclination and speed were calibrated using the Qualisys Pro Reflex 
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system and the Qualisys Track Manager software (Qualisys AB, Gothenburg, Sweden). 
The treadmill belt consisted of a non-slip rubber surface that allowed the skier to use 
his own poles (pole length: 90% of body height) with special carbide tips. The skier 
used a pair of Swenor skating roller skis with standard wheels (Swenor Roller skis, 
Troesken, Norway) and the Rottefella binding system (Rottefella AS, Klokkartstua, 
Norway), and the roller skis were pre-warmed before each test through 20-min of 
roller skiing on the treadmill. The roller skis were tested for rolling friction force (F f ) 
before the test period, and the friction coefficient ([i) was determined by dividing F f by 
the normal force (N) (|i = Ff • N" 1 ). This was performed in a towing test on three sub- 
jects (70, 80 and 90 kg) rolling passively at 3.9, 4.4 and 5.0 m/s for 5 min on a flat 
treadmill (0%) whilst connected to a strain gauge force transducer (S-type 9363, Revere 
Transducers Europe, Breda, The Netherlands). The measured [i was independent of 
speed and body mass, and the mean (i-value (0.0237) was incorporated in the work 
rate calculations. Gas exchange values were measured by open-circuit indirect calori- 
metry using an Oxycon Pro apparatus (Jaeger GmbH, Hoechberg, Germany). Before 
each measurement, the V0 2 and VC0 2 gas analyzers were calibrated using high-preci- 
sion gases (16.00 ± 0.04% 0 2 and 5.00 ± 0.1% C0 2 , Riessner-Gase GmbH & co, Lich- 
tenfels, Germany), the inspiratory flow meter was calibrated with a 3 L volume syringe 
(Hans Rudolph Inc., Kansas City, MO). Heart rate (HR) was measured with a heart 
rate monitor (Polar S610, Polar Electro OY, Kempele, Finland), using a 5-s interval for 
data storage. Blood lactate concentration (BLa) was measured on 5 [iL samples taken 
from the fingertip by a Lactate Pro LT-1710 t (ArkRay Inc, Kyoto, Japan). 

In a first test, the skier performed 5-min constant work rates at 0.05 and 0.12 
inclines in radians when treadmill roller skiing in the skating G3 technique. Gas 
exchange values were determined by the average of the last minute during each stage. 
The lactate threshold was defined at the metabolic power when blood lactate began to 
accumulate (OBLA) (defined as a concentration of 4 mmol/L, as calculated by a line- 
arly interpolated point out of the three measurement points of blood lactate concentra- 
tion at the incline of 0.05. Maximal metabolic power was tested at an incline of 0.05 in 
the G3 technique with a starting speed of 4.4 m/s. The speed was increased by 0.3 m/s 
every minute until exhaustion. V0 2 was measured continuously, and the average of the 
three highest 10-s consecutive measurements determined V0 2 max and used to calcu- 
late the maximal metabolic power. The test was considered to be a maximal effort if 
the following three criteria were met: 1) a plateau in V0 2 was obtained with increasing 
exercise intensity, 2) respiratory exchange ratio above 1.10, and 3) blood lactate con- 
centration exceeding 8 mmol/L. 

Thereafter, the velocity and the angle of incline on the treadmill were varied according 
to Figure 2 with gas exchange measured continuously. Immediately after finishing the 
protocol in Figure 2, the skier had a recovery period while skiing at 0.05 incline at 2.2 
m/s, inducing a work rate of 125 J/s and an aerobic power of approximately 900 J/s. 

Results 

The simulated data for steady state concentration of blood lactate during 5-min con- 
stant work rates show good agreement with the experimental results, with less than 0.5 
mmol/L disagreement (Figures 3 and 4). For the protocol showed in Figure 2, the 
simulated and experimental data for work rate and metabolic power as functions of 
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Figure 2 Velocity and incline as a function of time while treadmill roller skiing using the skating 

G3 technique. : The skier's velocity (v) in 10 2 m/s as a function of time (t) in seconds....:The treadmill's 

incline (a) in radians as a function of time (t) in seconds. 



time show general good agreement (Figure 5). However, we see that Q a (t) is delayed 
compared to the quasi steady state value Q a (t). Additionally, Q a (t) is delayed com- 
pared to the experimental data because the measuring apparatus of V0 2 has an inher- 
ent time lag (around 15 s) that is not modeled. Based on the same test (i.e., Figure 2), 
anaerobic power is simulated in Figure 6. Thereafter simulated and experimental lac- 
tate levels as functions of time were compared in Figure 7. There was a good agree- 
ment between the simulated and experimental results at these varying sub-maximal 
intensities, showing less than 0.5 mmol/L disagreements. During the recovery period 
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Figure 3 Steady state concentrations of blood lactate (C) in mmol/L as a function of the fraction of 
maximum aerobic power while treadmill roller skiing using the skating G3 technique at a 0.05 
incline. ■: Experimental data : c=£= i Ln (i + doQ °° t ^ Q ° ,J ) /;c - 
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after performing the protocol in Figure 2, the initial simulation model of the concen- 
trations of blood and muscle lactate did not fit well with experimental results. The 
experimental data showed a much slower reduction in blood lactate concentration 
compared to the simulated model. Thus, in a second trial we assumed a maximum 
rate of flux of lactate from blood to muscles. We achieved this mathematically when 
C m (t) - C b {t) — » Max[-0.1, C m (t) - C h (t)\ which accounts for a restriction on the lactate 
flux from blood to muscles. This means that the flux of lactate from the blood to the 




200 400 600 :jj 1000 1200 

time t 



Figure 5 Calculations of work rate (W) and metabolic powers (Q) as functions of time (t) in J/s for a 
skier while treadmill roller skiing in the skating G3 technique at the velocities and inclines shown 
in Figure 2. : Upper Q a (t), : Lower: Work rate (W), : Experimental data, :Q a (t). 



Moxnes and Sandbakk Theoretical Biology and Medical Modelling 2012, 9:7 
http://www.tbiomed.eom/content/9/1/7 



Page 11 of 14 



muscles has a least upper bound. The value of -0.1 was the value that fitted the experi- 
mental data best by visual inspection. Using this assumption, the simulated and experi- 
mental data were in much better agreement (see Figure 7). 

Conclusions 

The current study constructed mathematical models of lactate production and removal 
and compared these with experimental results from treadmill roller skiing. The main 
findings were as follows: 1) a mathematical two-compartment model of lactate produc- 
tion and removal could accurately predict blood lactate concentration during steady 
state and at varying exercise intensities and; 2) the understanding of lactate removal 
after high-intensity exercise during whole-body exercise requires further examination. 

In the current study, the simulated lactate production and removal fitted the experi- 
mental data well during steady state and at varying sub-maximal exercise intensities, 
and indicate an uncertainty of less than ±0.5 mmol/L using these mathematical mod- 
els. To the best of our knowledge, the current study is the first to compare a simula- 
tion model of blood lactate concentration during dynamic whole-body exercise against 
experimental data. Overall, we propose that the current simulation models provide 
useful insight into how blood lactate concentration during practical training and test- 
ing in dynamical situations should be interpreted. 

The initial simulation model of lactate removal after high exercise intensities 
deviated from the experimental findings. This can be related to a longer time delay 
between lactate concentrations in muscle and blood than simulated by the model [22] 
or a slower removal of blood lactate than the bi-exponential time function employed 
here [26-28]. In any case, when assuming a maximum rate of flux of lactate from 
blood to muscles did the simulation model fit experimental data satisfactorily. A ratio- 
nale for this asymmetry of the time scale for influx and outflux of lactate from the 
blood may be related to the transport of lactate across membranes accomplished by 




1750 



time t 

Figure 6 Simulated anaerobic power (Q) in J/s as a function of time (t) in seconds for a skier while 
treadmill roller skiing in the skating G3 technique at the velocities and inclines shown in Figure 2 
and a subsequent recovery period at 0.05 incline at 2.2 m/s. 
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Figure 7 Simulated blood and muscle lactate concentrations (C) in mmol/L as a function of time (t) 
in seconds for a skier while treadmill roller skiing in the skating G3 technique at the velocities and 
inclines shown in Figure 2 and a subsequent recovery period at 0.05 incline at 2.2 m/s. *: 

Experimental data, : Muscle compartment simulation, : Blood compartment simulation. 



monocarboxylate transport proteins (MTC). MCT are used for active transport of lac- 
tate and is driven by concentration gradients [9,13]. However, the physiological role of 
the two important isoforms MCT1 and MCT4 differ. MCT1 is mainly located in slow 
twitch muscle fibers and is regarded responsible for the influx of lactate to the high 
oxidative muscles, whereas MCT4 is mainly located in fast muscle fibers and is respon- 
sible for outflux of lactate from cells. The current data indicate that the amount of 
MTC4 is low in blood cells, which leads to a smaller outflux of lactate from the blood 
during recovery than the influx of lactate to the blood during high intensity exercise. 
Nevertheless, the understanding of lactate removal after high-intensity exercise during 
whole-body exercise requires further examination. 

Modeling in human biology is always a challenge since one is confronted with con- 
ceiving a simple, yet realistic representation of complex phenomena occurring at differ- 
ent levels (cells, organs, tissue). The parameters used in the current study will, in 
principle, be dependent on the exercise mode employed; fitness level of the individual 
tested and is only applicable for fixed concentrations of glycogen e.g., [42-44]. Thus, to 
construct valid simulation models of lactate concentration during dynamic exercise 
these need to be developed for each individual. 
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W, Work rate in J/s; v, Velocity of the treadmill in m/s; a, Angle of inclination in radians; m = 77.5 kg, Mass of the 
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Maximal aerobic power in J/s; QLT=0.93Qmax:, Aerobic power at the lactate threshold relative to the maximal aerobic 
power; Qa, Aerobic power in J/s; Qa, Steady state aerobic power in J/s; Qa, Virtual aerobic power in J/s; Q=a, Virtual 
steady state aerobic power in J/s; Qb = 80 J/s, Resting aerobic power; Q 0 (r 0 ) = Q b , Initial aerobic power; Qui = 111 J/ 
s, Aerobic power of unloaded movement of arms and legs; f, Cycle frequency of the skier in the G3 technique in 1/s; 
A(ct)=0.92(1+1.19Exp(-71.8a)), Function that describes the influence of incline on Qa; c 2 = 5.8 J/s, Parameter that 



Moxnes and Sandbakk Theoretical Biology and Medical Modelling 2012, 9:7 
http://www.tbiomed.eom/content/9/1/7 



Page 13 of 14 



describes the influence of work rate on Qa; r = 30 s, Time parameter in seconds quantifying the time before the 
aerobic power reaches steady state during sub-maximal work rates; C, Lactate concentration in the lactate pool in kg/ 
ml or in mmol/L; (C), Steady state lactate concentration in the lactate pool in kg/ ml or in mmol/L; C(f 0 ) = 0.045/cg/m 3 
= OSmmol/L, Initial concentration of lactate; R ai Rate of lactate appearance in kg/(sm 3 ) or in mmol/(sL); R dl Rate of 
lactate disappearance in kg/(sm 3 ) or in mmol/(sL); P max , Maximum rate of pyruvate appearance in kg/(sm 3 ) or in 
mmol/(sL); P, Rate of pyruvate appearance in kg/(sm 3 ) or in mMol/(sL); P, Steady state rate of_pyruvate appearance in 
kg/(sm 3 ) or in mmol/(sL); P, Virtual rate of pyruvate appearance in kg/(sm 3 ) or in mmol/(sL); (P), Virtual steady state 
pyruvate appearance in kg/(sm 3 ) or in mmol/(sL); r an = 10 s, Time parameter quantifying the time before the pyruvate 
appearance reaches steady state during sub-maximal work rates; a 0 = 0.9, Parameter that describe the rate of 
pyruvate disappearance due to oxidation in the mitochondria; /3 0 = 1/(0.6 x p 0 x Q max ), Parameter that describes the 
rate of pyruvate disappearance due to oxidation in the mitochondria; p 0 = 1 0~ 5 kg/{m 3 s)/{J/s), Parameter that describes 
the rate of lactate appearance as a function of aerobic power; D(Q 0 (r)), Function that describes the rate of lactate 
disappearance as a function of aerobic power; d 0 = 7.241 0~ 8 /C//s) 2 /s, Parameter that describes the rate of lactate 
disappearance as a function of aerobic power; \ = 0.95m 3 //cg, Parameter that describes the saturation of 
disappearance of lactate; QanG, Power caused by the use of ATP produced anaerobically from glycolysis or 
glycogenosis; A = 3 x 20J/{kgmmol/L), Parameter that scales the rate of change in lactate concentration and 
anaerobic power; V m = 10L, Lactate volume of the muscles; V b = 4L, Lactate volume of the blood; V = V b + V mi Total 
lactate volume; K ]a = 0.05/s x L, Parameter that scales the movement of lactate into and out of the blood 
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